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ABSTRACT 

Statistics are presented for the estimates of rain rate from an attenuating 
frequency radar. Included in the estimates arc the effects of fluctuations In 
return power and the rain- rate/reflectivity relationship, as well as errors 
introduced in the attempt to recover the unattenuated return power. In addi- 
tion to the Hitschfeld-Bordan correction, two alternative techniques are con- 
sidered. The performance of the radar is shown to be dependent on the method 
by which attenuation correction is made. 
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STATISTICS OF RAIN-RATE ESTIMATES FOR 


A SINGLE ATTENUATING RADAR 
I. INTRODUCTION 

Meteorological radars that require relatively high spatial resolution or 
sensitivity, but that are limited as to the antenna size or transmitted power, 
make desirable the use of higher frequencies. At X-band frequencies and higher, 
however, the signal may undergo significant attenuation by the intervening rain 
as it travels to and from the rain- scattering volume of interest. Therefore, 
the utility of such a system is in part dependent on the accuracy by which the 
nonattenuated power, as indicative of the rain rate in the scattering volume 
alone, can be constructed from the measured return power. 

To estimate the rain rate, R, for a single radar system, the reflectivity 
factor, Z t , is required. The latter quantity can be found directly from the 
radar equation only if the attenuation, k, is expressible as a function of either 
the reflectivity factor or some measurable quantity of the radar system. 

Because a number of such approximations for attenuation are possible, the 
question as to which rain- rate estimate should be employed reduces to the choice 
of the approximation that in some way optimizes radar performance. As will 
be shown, the accuracy of the rain- rate prediction is strongly dependent on the 
particular approximation adopted. 

The procedure given here is to form three different estimates of rain rate 
that involve the empirical (k, ZD, (Z t , R) relations and the measurable 
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quantities of the radar system. 1' renting as random variables the return power 
ironi each range bin, as well as certain parameters that relate (k, Z{) and 
(Z t , R), the mean and variance of each of the estimates are then computed as 
functions of range and number of indcitendent samples. Calibration errors 
and finite signal-to-noise ratios are also considered. 

II. RAIN-RATE ESTIMATES 

From the standard meteorological radar equation, including attenuation 
effects, the measured return power may be written as (Battan, 1 973] 

P(r) * CZ t (r) [l& 0 2 £ m*> (1) 

On using the relation k &Zj, then, 

hr) = {7,<r» [ia" J 11 J .’ J ]/r 3 (2) 

where k one-way attenuation in dB/kin 

P = measured radar return power averaged over X independent samples 
r = range in km 
C radar calibration constant 

For convenience, the auxiliary quantity, Z m , is introduced in terms of 
the defining relation: 

hr) • ( Z <r)/r 2 (3) 
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Or, on combining equations 1 and 2, 


Z,(r) ■ Z m <r» I0 0-2 a (4) 

The quantity, Z m , may be interpreted as a measured reflectivity factor 
in the sense that it can be found directly from the measured power and the known 
values of range and calibration constant. 

The reflectivity factor, Z ( , and rain rate, R, are usually given in the form 
of a power law relation that may be written as R - aZ^'. To estimate the rain 
rate, it is necessary to express Z t in terms of the measured quantity, Z m ; or, 
more simply, given the function, Z n<Z t ) (equation 1), it is required to find 

the solution of the inverse problem, Z, g _2 (Z m ). If such a solution exists, 
the rain rate then becomes F. - a [ g 1 (Z nl )] . 

To invert equation -t, it can Ik seen that, if all attenuation effects are ne- 
glected, a crude approximation io Z t would be 

°Z(r) = Z (r) (5’ 

in ' 

For any higher order approximation, the Z t that appears in the exponent 
of equation 4 is replaced by its preceding approximation. For first- and second- 
order iterations, this procedure gives 

'Z(r) = Z jr) 10"- 2 ° 6 ( 6) 
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An exact solution for Z t g - ' (Z m ) can be found by using a procedure 
similar to that given bv Hlischfeld and Iiordan (11)54). Taking the logarithm of 
both sides of equation 4, differentiating vvith respect to range, and, finally, 
using the substitution, u Z( , yields the first-order linear differential equation, 

dll , 

— + du Hr) + kd 1 0 

dr («) 

where 


then, 


k' = 0.2 q In 10 
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The condition, Z t (o) Z n (o), requires that c 1, and therefore, 







/•* <w)dw *' * 
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Approximating the integrals in the foregoing equations by summations over 
the relevant range bins :uid substituting the expressions for Z t of equations 
(>, 7, and 10 into the Z-R relation yield the following estimates of rain rate 
at the n ltl range bin: 


1 K aZ exp < Kh 5"! ' Z ) 

inn 1 L.~d i mi 


(in 
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J R = .1/ ' exp (Kb t / exp (-dk y* t //*. ) i 
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(12) 
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where 
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" "R = aZ b /() • kp T] c Z * ) h ^ 
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/Z** (s) ds ^ s Y* e Z*' 
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s = range resolution in km 


i=l i * n 

I 

i = n 
e = 1 j ^ i 

j 

i = i 
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The estimates of rain rati* given In equations 11 to 13 are oxact In the sense 

that statistical fluctuation ;tnd offset errors have been » \cludcd. II N lndc|>eudcnt 

th 

pulses are summed over the i lM range l)in, the average return power is a random 
variable given by (Marshall and Milschfeld, 1953; Stogrvn, 1975) 

or 


when the quantities, i’j and Z n) |, an the expectations of the return power and 
the measured reflectivity fact< . , respectively* rhe pro) rtlity density function 
of f is 


ptfl = 


N N 
IN - I |! 


c N, > ' 


(f >0) 


(14) 


If, in addition, calibration errors are introduced in which C C tru , II * E), 


then. 


/ = / . f /( I + It 

tm mi i 


(13) 
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Replacing Z ni j in equations 11 through 13 with '/ IT1 j results in the following 
rain- rate estimates: 
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where 


K * 0.: so In 10 

c * 1/(1 + fi) * C UU JC 


Quantities a and o In the above equations are chosen to be independent 
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random variables with variances of <* a and » 0 and means equal to their re- 
spective tr^e values. Quantities b and P are assumed to be fixed. 

Appendix A describes the behavior of the iterative approximations and their 
relationship to the Hitschfeld-Bordun solution. Appendix B provides a gener- 
alization of the estimates to the case of finite signal-to-noise ratios. 


*.fl. STATISTICS OF THE RAIN-RATE ESTIMATES 
The measures of performance for the rain-rate estimates are taken to be 
the mean and variance. Although it is sometimes necessary to use Monte Carlo 
techniques for generating the statistics, under certain conditions approximate 
analytic expressions can l>e found by direct integration. 

A. First-Order Estimate 

Assuming that the random variables (o, a, fj, f 2 , ... f n ) arc independent, 
the variance of the first-order estimate at the n th range bin, H^, may be 
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written as. 
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and 

ht.it « a , vjrt.r l * a 2 

m a 

p ( * 0.2 > c *' b In 1 0 

For large N, the expectati m with r« sped to t ( i l , 2 , ... . 
approximated by the first-order term of a saddle-point integration. Because 
(3 ^ usually near 1, an alternative procedure may be used that yields a some- 
what more accurate result. To evaluate the expectation with respect to 
fj ( i 1, 2, ... n-1), the term exp (i.' f.' ) (where’.. p.oZpjj is treated as 
a constant) is expanded about (3=1; for the integration over f n , exp (i;f ) is 

expanded alx>ut .dp 0. For the final integration over o, a normal distribution 
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is assumed with mean o m and variance r G . 

Neglecting terms involving powers of o greater than 2, the first and second 
moments of the first-order rain rate may be written 


i f 1 7 • ■ 1 £ k 
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where 
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D, = p, (fi - 1) (Q - In N) 

U 2 = P 2 (P - I )/N 

O = I + 'A + + l/N - 0.5772 1 57 
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H <q. p) 
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1(1 + 2q : n|/: P 
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* i 


H* = H. (q'. pi 


For the forenoinu unprtmed quantities, 6 Q 0; the primed quantities are given 
by the same expressions, but with 5 0 - 1. 

B. Second-Order Approximation 

The rain-rate estimate for the second-order approximation may be written 


in the form- 
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where. 


P 2 = 0.2 s fle* In It) 
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With respect to the independent random variables (fj, f 0 


f n ), the 


expectation in equation 25 may be written as 

Nn 


N 


; ff p 


[(N - I)’]" ” 

n I 

r^p, £ f . z i x f ex p <a Pj n s c/A * jx . i,x j • • dx „ 

*• « i * i j 

P'or large N, the saddle points occur at approximately x. 1 (i 1, 2, .... n). 
Evaluation of the integral results in the asymptotic approximation [Gradshteyn 
and Ryzhik, 1965; Felscn and Marcuvitz, 1973| 


II fc_ ( ) * exp 

i= i " 


«P [3^.2:*,^ ^ *' /m, ]/[l + l/l -N + 1/288N 2 ]” 


An identical procedure is used for equation 26. Taking the expectation with 
respect to o, letting u = p(o - o m ). and combining the previous results gives 


fc( R i * c J (o : + a 2 ) Z 

it j nt 


OB 

>Z l j e’ 0 * gfu) ( 1 + I/I2N * I/288N 2 )" 


(27) 


K(K ) % ca Z 

n tn m n 
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f e’ u *g'(ut du/^ir 

•/ ao 


(I + I/I2N + 1/288N ) (28) 


where 


g(u) = exp 


I 


g'(u) = ygdi) 


n 


The remaining integrals were computed by means of a Gaussian-Hermite 
quadrature. 

C. Hltschfcld-Bordan Estimate 

A problem arises in calculating the statistics of the Hitschfeld-Hordan 
estimate because the rain-rate estimate becomes imaginary when 

•»' £ *. < C > ' 

i= i 

To circumvent the problem, a Monte Carlo technique is used. As before, 
the estimate is characterized by its mean and variance. In addition, however, 
it is necessary to define a failure rate as the percentage of times the estimated 
rain rate is imaginary. 

IV. RESULTS AND DISCUSSION 

The results of the three rain-rate estimates for a radar operating at 8.75 
GHz are shown in the Figures. The Z t -R, k-Z t relations were derived from the 
tabulated values of Medhurst 11905] mid Stephens [1961] and are approximated 
by Z t = 307.1 k = 5.5 x 10"° Z ^ Sl . For reasons of economy, both the 

mean and variance are plotted on the same graphs as functions ol range for 
values of N of 100 and 1000. The "X” and "O" labeled points represent, re- 
spectively, the variance and mean of the rain-rate estimates (normalized to the 
true value) as computed by the Monte Carlo method. Where present, the solid 
and dashed curves represent the variance and mean as calculated from the 
approximate analytic expressions found in Section II. Unless otherwise noted, 
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an increase in N results in a decrease in variance and an increase in mean. 

The change in the latter quantity, however, is not significant. 

Comparison of analytic ;md Monte Carlo evaluations of the stall Hies are 
shown in Figures l and 1, and 2 anti 5 for the first- and second-order estimates, 
respectively. In addition to the qualitative insight that may be gained, the 
analytic results eliminate the spurious fluctuations that arise from the Monte 
Carlo generation at small values of variance. It ean be seen that the approx- 
imation for the second-order estimate for N = 100 is somewhat crude, implying 
that higher order terms are needed in the multiple-integral saddle-point tech- 
nique. 

The remaining graphs presented here were found by computing the sample 
mean and variance of 1000 simulations per range bin. To achieve adequate 
plotting resolution, the upper limit for the normalized variance was set at 2, 
which resulted in zero failure rates for the Ilitschfeld-Bordan estimates that 
were plotted. 

The most prominent features of the results follow: 

• In the absence of radar calibration errors and for low rain rates and 

small 'a , , the Hitschfeld-Bordan estimate is relatively unbiased- 

a consequence of the fact that the true reflectivity factor has been 
exactly expressed as a function of the measured reflectivity factor, Z 

• Even in the absence of calibration 2 rrors, the first- and second-order 
estimates are biased by an amount proportional to the total attenuation. 
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The use of higher order estimates tends to decrease Uiu offset error 


(Figures 1, 2, 4, and 5). 

• In general, the error variances are largest ;ind the sensitivity to 
calibration errors are most pronounced when a Hitschfeld-Dordan 
estimate is used (Figures 3 and (i). Because of the appearance of the 
ea and c^o factors in all the estimates, the effects of calibration 

errors are dependent on the statistics of a and o. Note that the 

o ^ 23 

variances of these factors are c 4 ° a , e . If c < 1 (i.e., the radar 
calibration constant is larger than the true value), a smaller error 
variance is expected than if the calibration constant were underpredicted 
(c > 1). This typo of behavior is indicated in Figures 7 through 12. 

• The rain-rate estimate for nonattenuating radar at the n 1 *’ range bin 
can be written as 

R „ 3 “ z !'„ r. 

The two characteristics that distinguish R n from the attenuated es- 
timates are: 

a. The calibration error enters into the estimate only through 
the multiplicative factor, c. 

b. The estimate is independent of the preceding range bins. 

A comparison between attenuating and nonattenuating radars can be obtained 
by noting that the mean and the variance of the latter are independent of range 
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(assuming largo signal-to-noiso ratios) and nearly equ; 1 to the values of the 
attenuated radar in the first range bln. The vertical deflection from this con- 
stant mean and variance indicates the additional errors that arise in using an 
attenuated frequency in conjunction with a particular reconstruction technique. 

V. CONCLUSIONS 

Three attenuation correction procedures have been considered that have 
led to the corresponding rain- rate estimates. For all estimates, an attempt 
has been made to include the effects of fluctuations in the averaged radar return 
power and the variability in the Z-R, k-Z relations. For several combinations 
of meteorological conditions and radar parameters, the mean and variance were 
computed and plotted as a function of range. 

The results indicate that, in addition to the Hitschfeld-Bordan procedure, 
alternative correction schemes exist that lead to rain-rate estimates with smaller 
error variances and a decreased sensitivity to positive offsets in the radar cal- 
ibration constant. On the other hand, corrections such as the first-order algo- 
rithm may lead to estimates that significiintly underpredict the true value of 
rain rate. 

Although it is impossible to choose the "best" correction technique without 
specifying the radar parameters and performance criteria, a few general 
characteristics of the various order estimates can be 3een. 

For a well-calibrated radar using relatively accurate k-Z, Z-R relations, 
higher order correction techniques may be employed. When the possibility of 
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significant errors exist, the lower order estimates must be used to avoid large 
error variances and large positive biases of the estimate. 

It is reasonable to assume that, given a radar design and performance 
criteria, there will be some order estimate (with k finite*) that optimizes 
the radar performance. For the criteria given, this estimate will represent 
the best compromise between the extreme cases of no attenuation correction 
and the attempt at exact correction. 
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APPENDIX A 


The n^ order approximation to Z t is given by 


n C r n • I d 

Z(r) = Z m (r)exp(K I Li r) 


dr) 


subject to the conditions that 


Z(r) = Z„, (r) 


Z(o) = Z (o) = Z, (o) 

m t 


(Al) 


Taking the logarithm of both sides of equation Al, differentiating with 

k -P 

respect to range, and letting u^ = Z yields the difference differential equa- 
tion, 


du (| /dr = - 0u, (f(r) + K/u n _ , ) 


(A2) 


where 


f(r) = d (In Z m )/dr = (1/Z m ) dZ m /dr 


From the definitions of Z(r) and Z^(r) and by induction, the following 
two conditions can be shown to hold for all values of r: 


n _ n • I _ 

Z(r) 5* Z(r) 
Z ( (r) > "Z(r) 
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I'o show that lltn n Z(r) - Z.(r), consider the series and an upper bound, 
n— • 

l,m 52 t /<r» - /m)=lim | N /ul / trl i t/tit / ml 
N N * oo m ' 


From the solution of the differential eqt alien (equation 10), it can In- seen 

that Z t (r) is bounded if K / Z |n (r)dr<l. Therefore, under this condition, 

Jo 

the series converges, implying that 


lim C L - " ' Z) * t) 

n -• 

or, equivalently, 

I mi * ll „/o ii , I = lllll ( // Z) ■ I 

n -*«• n -» °° 


In the limit of large n, the difference differential equation becomes 


du n /dr * -0ffr)u n - 0K (A3) 

subject to the condition, a Z(o) - Z m (o). i Jut this differential equation and 
initial condition are identical to those satisfied by Z ( (equation !)), so that 

lim Ztri = Z ( (r) 
n 
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APPENDIX 11 


The generalization of the estimates to the cuee of finite signal -to- noise 
ratios is straightforward. The power available at the receiver is mm the sum 
of N Independent samples of return power, plus receiver noise. Assuming that 
the noise power, p.,, Is statistically independent of the bach scattered power, 
p(r), but follows the same distribution, the power from the receiver can Ik- 
written as (Stogrvn, I !»7f>| 

P * (Ptr) + P N ) I <IU) 

As before, defining the measured reflectivity factor by 


/ 5- h r)r l /(' 

Ml 

where 

C ■ t, ru .d ' Li 

then, 

Z m = (Plri ♦ P N >fr 2 /C |ru ,<l ♦ l » 

or at the i**' range bin, 

Z =7 f (I + I* IP (r))/( I + I l (B2> 

Comparison of equation 112 with equation 1 > of the main report shows that 
the effect of finite signal to noise is to introduce an additional term (1 + p_/p(r)) 
into the definition of the measured reflectivity factor. Proceeding as Indore, 
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Z t Is expressed as a function of Z m by various order approximations that. 
In turn, lead to three estimates of ruin rate analogous to those of equations 




16 through 18, 

For example, the first-order rain rate becomes 
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(B3) 
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Figure 9. Statistics of the Rain-Rate Estimate for Hltschfeld-Bordan 

Attenuation Correction (R = 25 mm/hr, ® a 3 a/20, 0 o = °/10, c - 1.25). 

Fifmre 10. Statistics of the Rain-Rate Estimate for First-Order Attenuation 
Correction (R - 25 mm/hr, o a = a/20, <* 0 * o/ 10, c = 2/3). 


Figure 11. Statistics of the Rain-Rate E timatc for Second-( 'rder Attenuation 
Correction (R = 25 mm/hr, ^ a = u/20, <> 0 = ®/10, c = 2/3). 

Figure 12. Statistics of the Rain-Rate Estimate for Hitschfeld-Bordan 

Attenuation Correction (R 25 mm/hr, <r a = a/20, <* Q o/10, c = 2/3). 
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Figure 3. Statistics of the Rain-Rate Estimate for Hitschfeld-Bordan Attenuation Correction (R = 25 mm hr 

<% = a/10, o a = a/5, c = 1). 
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Figure 6. Statistics of the Rain-Rate Estimate tor Hitschfeld-Bordai Att» , 

«R - 40 mm/hr, , . a , 20. = itV = T °" C ° rreCti ° n 
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I igure 7. Statistics of the Rain-Rate Estimate for Tirst-Order Attenuation Correction 

(R = 25 mm/hr, ? a = a, 20, <% = 10, c = 1.25). 
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Figure 8. Statistics of the Rain-Ilate Estimate for Second-Order Attenuation Correction 

(R = 25 mm/hr, cr # = a/20, ~ a = a/10, c = 1.25). 
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Figure 9. Statistics of the Ram-Rate Estimate for Hitschfeld-Bordan Attenuation Correction 

(R = 25 mm/hr , - a/20, % = x/10, c = 1.25). 
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Figure 10. Statistics of the Rain-Rate Estimate for First-Order Attenuation Correction 

(R = 25 mm/hr, <r m = a/20, *'« = l / 10, c = 2/3). 
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Figure 11. Statistics of the Rain-Rate Estimate far Second-Order Attenuation Correction 

(R= 25 inm/hr, : t = a. 20, a = i/’lO, c = 2/3). 
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.-urc 12. Statistics of the Rain-Rate Estimate for Hitschfeld-Bordan Attenuation Correction 

(R = 25 mm/hr, c = a/20, c a = j / 10. c - 2,3). 
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